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ABSTRACT 

From optical photometry we show that SDSSJ121258. 25-012310.1 is a new eclipsing, post common-envelope binary with an orbital 
period of 8.06 hours and an eclipse length of 23 minutes. We observed the object over 1 1 nights in different bands and determined the 
ephemeris of the eclipse to HJDn,id = 2454104.7086(2) + 0.3358706(5) x E, where numbers in parenthesis indicate the uncertainties 
in the last digit. The depth of the eclipse is 2.85 ± 0.17 mag in the V band, 1.82 + 0.08 mag in the R band and 0.52 + 0.02 mag in the 
/ band. From spectroscopic observations we measured the semi-amplitude of the radial velocity ^"2 = 181 ± 3 km/s for the secondary 
star. The stellar and binary parameters of the system were constrained from a) fitting the SDSS composite spectrum of the binary, b) 
using a A'-band luminosty-mass relation for the secondary star, and c) from detailed analyses of the eclipse light curve. The white 
dwarf has an effective temperature of 17700 ± 300 K, and its surface gravity is log^ = 7.53 + 0.2. We estimate that the spectral 
type of the red dwarf is M4 + 1 and the distance to the system is 230 ± 20 parsec. The mass of the secondary star is estimated to be 
in the range M^^ = 0.26 - 0.29 Mq , while the mass of the white dwarf is most likely M„d = 0.46 - 0.48 Mq. From an empirical 
mass-radius relation we estimate the radius of the red dwarf to be in the range 0.28 - 0.31 Rq , whereas we get R^^ = 0.016 - 0.018 
Rq from a theoretical mass-radius realation. Finally we discuss the spectral energy distribution and the likely evolutionary state of 
SDSS 1212-0123. 
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1. Introduction 

Most stars are found in binary or multiple star systems and a 
large fraction of binaries will interact at some point in their lives. 
Interaction depends on initial separation, relative masses of the 
I components, and their evolutionary state. When the initial or- 
> bital period is less than ~ 5 days, interaction will take place 
while both stars are on the main sequence, giving rise to con- 
tact binaries and eventualy mergers. For wider orbits, interaction 
can happen only when the more massive star reaches the RGB 
or the AGB. If the mass transfer rate is high enough a gaseous 
envelope may surround the binary, entering thus a common en- 
velope (CE) phase. The stars spiral towards each other, friction 
within the CE will lead to a shrinkage of the binary separation, 
and angular momentum and energy are extracted from the or- 
bit, expelling the CE. Once the envelope is expelled, the remain- 
ing system consists of a main sequence star and a remnant of 
the more massive star, e.g. a white dwarf, perhaps surrounded 
by the ejected material that can be ionized f orming a planetary 
nebula (* Paczviiski|[T976t llben & Liviolll993h . Further evolution 
of the system is driven by angular momentum loss due to mag- 
netic braking and/or gravitational radiation, which will bring the 
system into a semidetached configuration, becoming a cataclis- 
mic variable (CV). Post common envelope binaries (PCEBs) are 
thought to be precursors of other interesting objects such as C Vs, 



* This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile. 



millisecond pulsars, low-mass X-ray binaries or double degen- 
erate white dwarfs. 

Though theories for the CE phase exist it is still poorly 
understood. The efficiency of mass ejection in the CE with 
respect to the masses of the components, their evolutionary 
state and the orbital separat ion is uncertain jPolitano & Weileii 
2007] iTaam & Ricked l2006|) . White dwarf plus main sequence 
(WDMS) binaries are a perfect test laboratory for study- 
ing the current population synthesis models because they are 
numerous, their stellar components are well understood in 
terms of their single evolution and they are not accreting, 
which would increase the complexity. The number of PCEBs 
with well-defin ed parameters is still small, although it has 
increased since ISchreiber & Gansickd (|2003l) from 30 to al- 



most 50_system^(Gansicke_eLal^^ Morales-Rueda et ajj 
2005 1; 'Shimanskv et al.j_i 2006^. Van den Besselaar et al. ,2007 



Aung weroiwit et aL 2007nT appert et al. 2007; Schrei ber et alj 



'20081; IRebassa-Mansergas et al.ll2008t ISteinfadt et al.ll2008 



In this paper we report the discovery of a new eclips- 
ing PCEB. In our ongoing search for PCEBs among 
white-dwarf/main-sequen c e biti ar ies jSchreiber et all 120081 ; 
iRebassa-Mans ergas et al.' 12007 1 |2008 |), SD SSJ121258.25- 
012310.1 iAdelman-McCarthy^talJ l2008[) (henceforth 
SDSS1212-0123) was included in our target list for pho- 
tometric monitoring of candidate objects. The serendipituous 
discovery of a binary eclipse from time-resolved differential 
photometry triggered a photometric and spectroscopic follow- 



2 



A. Nebot Gomez-Moran et al.: Post common envelope binaries from SDSS. IV: SDSS1212-0123, a new eclipsing system 



up. Only seven eclipsing binaries containing a white dwarf 
and a low mass main sequence star were known until 2007. 
Since then another three eclipsing syste ms have been published 
jSteinfadt et al.ll2008l:lD"rake etaLl2008l). and a furthe r three sys- 
tems have been discovered by us ( Pvrzas et al.ll200^ . Eclipsing 
binaries are of great interest since they offer the possibility of 
deriving fundamental properties of stars with a high accuracy. 
SDSS1212-012 3 was firstly listed as a quasar candidate by 
Richards et alj 
;t all 



Silvestri et ; 



and later classified as a DA + dMe by 
. It contains a relatively hot white dwarf 
(from now on primary) and an active mid-type dM star (from 
now on secondary). 

In this paper we summarize our current knowledge about this 
source from own observations and archival work. It is organized 
as follows. In Sect. |2] we describe the observations and reduc- 
tions. In Sect.|3]we present the results, we study the evolution of 
the system in Sect.|4]and conclude in Sect.|5] 

2. Observations and reductions 

2.1. IAC80 and AIP70 photometry 

Optical photometric observations were obtained using two dif- 
ferent telescopes over 11 nights. The 80 cm telescope IAC80 
in Observatorio del Teide, Spain, was equipped with the stan- 
dard CCD camera and the 70 cm telescope of the Astrophysical 
Institute Potsdam at Babelsberg was used with a cryogenically 
cooled 1x1 k TEK-CCD. A log of observations is presented in 
Table [T] A field of ~ 3 arc minutes was read with the IAC80 
CCD camera, and we used a binning factor of 2 in both spa- 
tial directions (scale of 0.6"), while we used a binning fac- 
tor of 3 for the 70 cm telescope (scale of 1.41"), in order 
to decrease the readout time and improve the signal to noise. 
Reduction was performed using standard packages in IRAlQ 
and MIDAS. Differential magnitudes were obtained with respect 
to the comparison star SDSS J121302.39-012343.5 (see FiglB, 
with magnitudes ugriz^ll. 40, 16.00, 15.51, 15.36, 15.30. SDSS 
magnitudes were transformed into Johnson's using equations 
taken from the Sloan pagefl Neglecting the color term, we cal- 
culated absolute magnitudes of SDSS1212-0123. The estimated 
error of the absolute calibration is 0.05 mag. 



Table 1. Log of photometric observations for SDSS1212-0123. 



Date 


Tel 


Filter 


tint 


A'obs 


0ini 


0fiii 


04 Jan 2007 


IAC80 


1 


180 


74 


0.805 


1.278 


26 Jan 2007 


AIP70 


V 


180 


19 


0.041 


0.156 


13 Feb 2007 


IAC80 


V 


70 


191 


0.607 


1.202 


14 Feb 2007 


IAC80 


V 


70 


221 


0.525 


1.292 


12 Mar 2007 


AIP70 


I 


120 


181 


0.031 


0.946 


13 Mar 2007 


AIP70 


1 


120 


153 


0.377 


1.343 


14 Mar 2007 


AIP70 


I 


120 


49 


0.786 


1.040 


15 Mar 2007 


AIP70 


R 


120 


73 


0.683 


1.023 


26 Mar 2007 


AIP70 


I 


120 


45 


0.994 


1.179 


21 May 2007 


AIP70 


1 


90 


59 


0.907 


1.097 


06 May 2008 


IAC80 


R 


120 


32 


0.913 


1.049 



' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under contract with the National 
Science Foundation, http://iraf.noao.edu 

^ http://www.sdss.org/dr7/algorithms/sdssUBVRITransform.html 
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Fig. 1. SDSS image of SDSS1212-0123 (in the cross-hair) and 
the comparison star (RA = 12:13:02.39, DEC = -01:23:43.5). 



2.2. Spectroscopy 

Spectroscopic follow up observations were obtained during the 
period 16-19 May 2007 with the LDSS3 imaging spectropho- 
tograph at the Magellan Clay telescope. Ten spectra were taken 
for SDSS1212-0123. Exposure times varied from 300 to 600 
seconds. Seeing and transparency were highly variable. The 
VPH_Red grism and an OG590 blocking filter were used. The 
detector was a STA 4kx4k pixel CCD with two read out ampli- 
fiers. We used a slit width of 0.75 arcsec, that together with the 
spectral resolution R = 1810, gave a coverage of 5800 - 9980 
A at a reciprocal dispersion of 1.2 Apix"'. Four of the spec- 
tra taken at quadrature were obtained through a narrow slit of 
0.5 arcsec resulting in a FWHM spectral resolution of 4.8 A, 
with the purpose of measuring the radial velocity amplitude 
with a higher accuracy. Flat-field images were taken at the 
position of the target to allow effective fringe removal in the 
red part of the spectra. The spectral images were reduced us- 
ing STARLINK packages FIGARO and KAPPA, and the spec- 
tra were optima lly ext racted (Horne 1986) using the PAMELA 
package (Ma rsh|[l989l) . Wavelength calibration was done using 
sky lines. Wavelengths of good sky lines were obtained from 
the atlas of lOsterbrock et al.l (|1996[ [1997). A fifth-order poly- 
nomial was fitted to 36 sky lines. Spectra were flux calibrated 
and corrected for telluric lines using spectra of the standard star 
LTT3218 taken during the same observing run. 



3. Results 

3. 1 . The light curve 

The optical light curve of SDSS1212-0123 displays a total 
eclipse of the primary with length of approximately 23 min- 
utes. The depth of the eclipse is 0.52 + 0.02 mag in the / band, 
1.82 + 0.08 mag in the R band and 2.85 + 0.17 mag in the V 
band (see Fig.|2]l. Eclipse magnitudes are m/ = 16™56 + 0^02, 
niR = 18™58 + 0'."08 and my = 19'J'68 + 0™17. The much deeper 
eclipse in the V band is due to the fact that the primary emits 
most of the light in the blue, while the secondary dominates in 
the / band. Photometric variability outside of the eclipse, e.g. 
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Fig. 2. Optical photometry from the IAC80 telescope in the V, 
R and / band (from bottom to top) phase folded over the orbital 
period. Note the different scales for each panel. 



from an irradiated secondary or from ellipsoidal modulation of 
the secondary, was found to be less than (FOl. 

At the given time resolution of our photometry, the WD 
ingress and egress phases are not resolved. Five eclipses were 
completely covered and the eclipse length was determined in 
these light curves measuring their full width at half maximum 
of the flux level. The weighted mean of those five measurements 
gives an eclipse length of 23 + 1 min. 



3.2. Ephemeris 

In addition to the five eclipses which were covered completely 
one further eclipse was covered partially. Using the measured 
eclipse length from the previous section we thus determined six 
eclipse epochs (Table |2]). The eclipses of March 12, 14 and 26, 
respectively, were not covered due to bad weather conditions. 
Using a phase-dispersion minimization technique a tentative pe- 
riod was determined, foib = 0.3359 + 0.0006 hours, which was 
sufficiently accurate to connect all follow-up observations with- 
out a cycle count alias. 

We then used the six mid eclipse epochs to calculate a linear 
ephemeris by fitting a line to the cycle number and eclipse epoch: 



HJD^d = 2454104.7086(2) + 0.3358706(5) x E, 



(1) 



where numbers in parenthesis indicate the Icr uncertainty in the 
last digits. The observed minus calculated values are tabulated 
in Table |2] 



Table 2. Date, times of mid eclipses, cycle number obtained 
from the photometric observations and residuals from the linear 
ephemeris. 



Date 


HID (Mid-eclipse) 


Cycle 


- C (s) 


04 Jan 2007'- 


2454104.7085(21) 





-9.9 


13 Feb 2007" 


2454145.6847(8) 


122 


-Il.O 


14 Feb 2007" 


2454146.6922(8) 


125 


-19.0 


13 Mar 2007" 


2454173.5621(10) 


205 


-0.3 


21 May 2007" 


2454242.4163(21) 


410 


66.9 


06 May 2008" 


2454593.4000(14) 


1455 


-25.9 



' Eclipse completely covered. 
'' Eclipse partially covered. 



3.3. Stellar parameters 

3.3.1 . Decomposition of tine SDSS spectrum 

We determined the stellar parameters of SDSS1212-0123 from 
the SDSS spectrum following the procedure described in 
iRebassa-Mansergas et al.l (|2007|) 

In a first step the best match of the SDSS composite spec- 
trum is determined with an optimization strategy on a grid of ob- 
served white dwarf and M-dwarf template spectra created from 
the SDSS DR6 database. The main result of this first step is 
the determination of the spectral typ e of the secondary. Using 
the sp ectral type-radius relation from iRebassa-Mansergas et al.l 
( |2007|) and the apparent magnitude of the scaled template re- 
sults in a first distance estimate t/sec- After subtracting the best- 
fitting M-star template, white-dwarf parameters are determined 
via minimization in a log g - T^ff grid of model atmospheres 
dKoester et alj|2005l) . Since this analysis step is performed on 
spectra normalized to a continuum intensity, the results are bi- 
valued yielding a 'hot' and a 'cold' solution (see Fig.|4li. The de- 
generacy can typically be broken by an additional fit to the over- 
all spectrum (continuum plus lines in the wavelength range 3850 
- 7150 A). In the present case of SDSS1212-0123 the GALEX 
detection (see below) provides an additional constraint exclud- 
ing the 'cold' solution. The results of the spectral decomposition 
and the white dwarf fit for SDSS1212-0123 are shown in Fig.|4] 
Mass and radius of the white dwarf are calculated with the 
best-fitting log g - T^ff combination using updated versions of 
the tables by Bergeron et al. ( 1995). The flux scaling factor to- 
gether with the derived radius of the white dwarf results in a 
second distance estimate of the binary, d„d. 

The spectral type of the secondary was determined to be 
M4 ± 1 implying a distance c/sec - 320 ± 95 pc, mass range 
of the secondary Msec = 0.25 5 - 0.380 Mq and radius range 



/?sec = 0.258 - 0.391 Rq, using lRebassa-Mansergas et all (l2007h 
spectral type-mass and spectral type-radius empirical relations 
respectively. The derived temperature and log g of the primary 
were found to be only weakly dependent on the chosen spec- 
tral type and spectral template of the secondary, because we 
use Hj6 - He for the white dwarf line fit, where the secondary 
star contribution is small. It is also weakly dependendt on the 
accuracy of the spectral flux calibration and also the small ra- 
dial velocity line displacements. The best fit was found for 
Teff = 17700 + 300 K and log g = 7.53 ± 0.05 (implying a white 
dwarf mass Mwd = 0.39 + 0.02 Mq, and R„d = 0.018 + 0.001 
Rq). The obtained values are in agreement with those published 
by Silvestri et al. (2006). However, one should be aware of the 
fact that all the quoted errors are purely statistical. The true un- 
certainty of the white dwarf spectral parameters is clearly higher 
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than suggested by the derived numbers. We estimate the system- 
atic uncertainty of our log g determination to be on the order of 
0.2 dex, which results in rather wide ranges of possible values 
for the mass and the radius of the primary, i.e. M^d = 0.33 - 
0.48 Mo and R„d = 0.015 - 0.021 Rq. 

The derived distance to the white dwarf is t/^d = 226 ± 8 pc 
(assuming the statistical error only). The two distance estimates 
differ, dsec being longer than d„d, but in agreement within the 
erros. lRebassa-Mansergas et alj ( l2007h found a similar trend for 
101 WDMS binaries in their study. They argue that such differ- 
ence could be due to stellar activity of the secondary star, and 
that the spectral type determined from the optical SDSS spec- 
trum is too early for the mass of the secondary star, which would 
lead to a larger radius and consequently a larger distance to the 
system. Since the secondary in SDSS1212-0123 was found to 
be active too, we regard the distance estimate for the white dwarf 
being more reliable. Taking into account systematic errors we 
obtain c/^d = 230 + 20 pc as the distance to the system. 

3.3.2. Constraining tine secondary mass using 2I\/IASS 

In the previous section we derived the mass and the ra- 
dius of the seco ndary star using empirical relations from 
iRebassa-M anser gas et alj ( |2007 ) and obtained Msec = 0.255 - 
0.380 Mo and rI^^ = 0.258 - 0.391 Rp, respectively. However , 
as clearly shown in Fig. 7 of iRebassa-Mansergas et all (|2007[) . 
the masses and radii derived from observations largely scat- 
ter around the empirical relations. In addition, according to 
iReb assa-Manserg as et al. ( 2007) increased activity of the rapidly 
rotating secondary stars in close binaries can cause the stars to 
appear as earlier spectral types when compared to non-active 
stars of the same mass. To sum up, the secondary masses de- 
rived from empirical relations can obviously only considered to 
reasonable but rough estimates. 

An alternative method to determine the mass of secondary 
star is to use luminosity-spectral type relations. To that end, 
we expl ored the Tw o Micron All Sky Survey Point Source 
Catalog dCutri et al.l2003b . finding magnitudes J = 14.90 + 0.03, 
H = 14.39 + 0.05 and K, = 13.96 + 0.05 for SDSS1212-0123. 
Subtracting the extrapolated contribution of the primary star 
(logg - 7.5 and c/wd = 230 pc) yields infra-red colors of 
J-H = 0.51 +0.06, H-K, = 0.43 ±0.07, respectively. U s ing th e 
empirical mass-luminosity relation from iDelfosse et all (|2000|) . 
we derive the mass of the secondary st ar to be 0.26 + 0.03. Using 
again the mass-radius relation from IRebassa-Mansergas et al.l 
( 1200 7) this implies a spectral type M5, i.e. later by one stectral 
type than estimated from the deconvolution of the SDSS spec- 
trum. This supports the idea of activity significantly affecting 
the determination the secondary star spectral types and the cor- 
responding distances. 

3.3.3. Radial velocity 

In each of our observed spectra we measured the radial veloci- 
ties of the Nal absorption doublet (8 183.27 A, 8 194.81 A), which 
originates from the secondary star. A double Gaussian with a 
fixed separation of 1 1 .54 A was fitted to the line profiles using 
the FIT/TABLE command provided by ESO/MIDAS. 

Ha was deconvolved into an absorption and an emission 
line component using two Gaussians. While the emission line 
showed pronounced wavelength shifts, the centroids of the ab- 
sorption lines thus measured did not constrain the curve of the 
white dwarf significantly. 
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Fig.3. Two component fit to SDSS1212-0123. The top panel 
shows the WDMS spectrum (black line) and the white dwarf and 
the M4 M-dwarf templates (dotted lines), while the lower panel 
shows the residuals to the fit. 
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Fig. 4. Spectral model fit to the white dwarf component of 
SDSS 1212-0123, obtained after subtracting the best-fit M-dwarf 
template from its SDSS spectra. Top left panel: best fit (black 
lines) to the normalized H/3 to He line profiles (gray lines, top 
to bottom). Top right panel: 1, 2 and 3cr contour plots in the 
^eif - logg plane. The black contours refer to the best line profile 
fit, the red contours to the fit of the whole spectrum. The dashed 
line indicates where the maxima of the H/3 equivalent width oc- 
curs in the Tgff-logg plane, dividing it into two different so- 
lutions, a cold and a hot one. The best-fit parameters of the hot 
and the cold normalized line profile solutions and of the fit to the 
3850 - 7 150 A range are indicated by the black and the red dots, 
respectively. Bottom panel: the white dwarf spectrum and asso- 
ciated flux errors (gray lines) along with the best-fit white dwarf 
model (black lines) to the 3850 - 7 150 A wavelength range (top) 
and the residuals of the fit (gray line, bottom). 

Assuming a circular orbit a sine-function was fitted to the 
measured radial velocity curves to obtain the radial velocity 
semi-amplitude K2 of the secondary star: 



Vr -72 + K2 sin 



(2) 
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Fig. 5. Radial velocities measured from the Nal dou- 
blet 8183, 81 94 A originating from the secondary star of 
SDSS 1212-0123 folded over the orbital period obtained from 
the photometry. Sine fit and residuals (lower panel) are shown. 

The orbital period P and the epoch of mid eclipse to were de- 
termined photometrically and were kept fixed for the radial ve- 
locity fit. For the Nal doublet we find the systemic velocity 
72 = 17 ± 3kms ' and K2 = 181 + 3kms"', while we find for 
the HcK line 72 = 21 + 2kms"'and /:2 = 161 + 3kms"'. The fit 
to the Nal lines is shown if Fig. |5] together with the residuals. 

The semi-amplitudes of the two radial velocity curves are 
different and these differences seem to be significant. The semi- 
amplitude derived from Ha is lower, indicating that its emission 
is displaced towards the inner hemisphere of the secondary star 
with respect to the Nal doublet. As neither of the two line fea- 
tures shows significant photometric variability, which would in- 
dicate a biased origin of one of the line species (e.g. towards 
the non-irradiated side of the secondary), we exclude irradiation 
as the explanation for the observed difference in K2- A detailed 
comparison of radial velocities deriv ed from the Nal doublet 
and H g lines has been performed by iRebassa-Mansergas et alj 
(1200 7'). They find that both velocities often significantly dif- 
fer but that there seems to be no systematic shift of Ha ra- 
dial velocities towards smaller v alues. As discussed in detail in 
iRebassa-Manserpas et alj ( l2008h . this is probably explained by 
the Ha emission being related to activity and not uniformely dis- 
tributed over the surface of the secondary. Kafka et al.l (l2005h 
studied in detail the origin of different line species, however 
SDSS 12 12-012312 12 shows no evidence of accretion nor ir- 
radiation. We therefore assume that in SDSS1212-0123 the Nal 
doublet much better traces the center of mass of the secondary 
and we use its semi-amplitude for the mass estimate. 

We write the mass function of the binary assuming a circular 
orbit in the form 



2;7rG sin^ / 



PKl 



-Mwri - 1 



(3) 



and derive an upper limit for Msec for a given white-dwarf mass 
M„d assuming ; = 90° (see bottom panel of Fig. |6l). 

Using the e mpirical mass-radius rela tion for main sequence 
stars derived by Bayle ss & Orosz (2006) we estimate the radius 
of the secondary (middle panel of Fig. |6]l. The top panel of the 
same figure illustrates the maximum possible eclipse length (/ - 
90°, black line) for the given stellar radius, the orbital period P 
and the orbital separation a according to 




fp.rl - 



(4) 



Fig. 6. Solution of the mass function for K2 - 181+3 km/s as 
a function of M„d (bott om panel); radius f rom th e mass-radius 
empirical relation from iBavless & Orosj ( 1200 6') (mid panel); 
eclipse duration (top panel), assuming circular orbit and an in- 
clination angle of 90°. Vertical lines indicate the mass of the 
primary star, M„d - (0.33 - 0.48) M©, as determined from 
the deconvolution of the SDSS spectrum. The eclipse length, 
feci = 23 + 1 min, is marked with horizontal lines, and the in- 
tersection in shaded in grey. The erros in K2 are shown with the 
small dashed lines in the three panels. In the top panel the solu- 
tion for an inclination of 75° is plotted with a blue line to show 
it's influence. See text for a more detailed description. 



The measured values of the eclipse length and the range of 
the white-dwarf mass from Sec. 13. 3.11 are shown in the figure 
with horizontal and vertical lines, respectively, their intersection 
is shaded in grey in the top panel. It is also plotted the solution 
for / - 75° for comparison. From the eclipse length the range of 
possible values for the mass of the WD is M„d - 0.46 - 0.52, 
and for the dM M^ec = 0.21 - 0.32 Mq, and R,^^ = 0.23 - 0.34 



3.3.4. Light curve modeling 

A determination of most of the physical parameters of an eclips- 
ing system can be achieved by fitting model light curves to the 
actual data. We made use of a newly developed light curve fit- 
ting code, written by T.R. Marsh, for the general case of binaries 
containing a white dwarf. 

The code is described in detail in Pvrzas et al.l ( l2008l) . 
Briefly, a model light curve is computed based on user-supplied 
initial system parameters. These are the two radii, scaled by the 
binary separation, R„d/a and /?sec/a, the orbital inclination, /, 
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the unirradiated stellar temperatures of the white dwarf and the 
secondary star Teff wD and Teff sec respectively, the mass ratio 
q - Msec/A^wd and fo the time of mid-eclipse of the white dwarf. 

Starting from this parameter set, the model light curve is 
then fitted to the data using Levenberg-Marquardt minimisation. 
Every parameter can either be allowed to vary or remain fixed, 
during the fitting process. 

Our approach for modeling the / band photometry of 
SDSS1212-0123 was the following. A large and dense grid of 
points in the M„d - Msec plane was first calculated, generously 
bracketing the estimates for the mass of the two components (see 
Sec. l3.3.3l l. Each point defines a mass ratio q, and through Porb, a 
binary separation a. Furthermore, from the mass function equa- 
tion (Eq. |3]l, using the value of K2 (derived in Sec. l3.3.3l l and 
Poib, one can calculate the inclination angle /. Points for which 
(formally) sin / > 1 were discarded from the grid, for all other 
points the corresponding light curve model was computed, lead- 
ing to the computation of some 9000 models. 

As an initial estimate for the radii of the binary compo- 
nents, we adopted values from the theoretical M - R relations 
of Ber geron et al (1995) for the white dwarf, and ^ Baraffe et al.l 
( Il998l) - the 5 Gyr model - for the secondary. Regarding the two 
temperatures, Teff.wo and Teff.sec , the value from our spectral 
decomposition (Sec. l3.3.l1 i was used for the white dwarf, w hile 
the Sp(2) - T relation from lRebassa-Mansergas et al.l(l2007l) . to- 
gether with our result for the spectral type of the secondary, were 
used to obtain an initial value for Teff^sec- 

For the fitting process q, i, R^^ and Teff wD were fixed, 
leaving only /?sec, 7^eft,sec and fo free to vary. 7?„d was fixed 
mainly because of the poor temporal resolution of our data 
set, which does not resolve the white dwarf ingress and egress. 
Consequently, if allowed to vary, the white dwarf radius would 
only be loosely constrained and it would introduce large uncer- 
tainties in the determination of 7? sec- Teif wD was also fixed, be- 
cause allowing both temperatures to vary simultaneously would 
lead to a degenerate situation, as they are strongly correlated. 
Our spectral decomposition results are sufficiently accurate, so 
as to allow us to fix Teff.wD without affecting the fitting result. 
The parameter to on the other hand, was left free during the fit- 
ting, to account for the O - C errors in the mid-eclipse times, 
which in some cases were significant (see Table |2] again). 

The results of the light curve fitting process were analyzed 
as follows. We first applied a cut in the quality of the fits. This 
was done by selecting the minimum value of all fits and then 
culling all model fits at > Icr above the best fit. Afterwards, we 
selected from the remaining, equally good light curve fits, those 
which where physically plausible. We defined a 6R parameter, 
as 6R - (Rfit - Rth) /Rth, i-e- how much the fitted radius value 
deviates from the theoretical radius value, obtained from a M-R 
relation, for a given model. Thus, we selected only those models 
that had 6R < 0.15, to allow for an oversized secondary. 

The results are illustrated in Figure|2] Black dots designate 
those light curve fits making the Icr cut, red dots those that 
satisfy both the Icr and 6R - 15% cuts. The resulting ranges 
in white dwarf masses and secondary star masses (indicated 
with dashed, vertical, red lines) are M^d - 0.46 - 0.6 and 
Msec - 0.23 - 0.4 Mq, respectively, corresponding to a white 
dwarf radius of ^wd = 0.013 - 0.016 R© and a secondary radius 
of /?sec = 0.27 - 0.41 Rq. The range for the inclination angle is 
i = 82° - 90°. Also indicated, with dotted, horizontal, gray lines 
are the radii of M-dwarfs with spectral types Sp(2) = M3 - M5 
in steps of 0.5, based on the spec tral type-mass relation given by 
iRebassa-Mansergas et alj (l2007h . 
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Fig. 7. Light curve model fitting results for SDSS 1212-0123. 
Left panel: Mwd and Msec values corresponding to fits with 
values within Icr of the minimum value (black points) and, si- 
multaneously, with 6R < 0.05 (red points). Right panel: the 
same, only in the q - / plane. Also depicted in the left panel 
are curves corresponding to the mass function (solid black lines, 
/ - 90° and 75°) which (by definition) bracket the possible so- 
lutions, Sp(2) - M relations (dotted, horizontal, gray lines) and 
the range of possible (Mwd, Msec) values (dashed, horizontal and 
vertical, red lines) 

FigUl shows one example of the light curve fits within the 
components masses range for the model parameters: M„d - 0.49 
Mq, Msec - 0.26 Mq and / - 89.2°. The detailed models do not 
predict any variation in the light curve caused by irradiation of 
the secondary star by the white dwarf. The predicted variations 
due to ellipsoidal modulation are expected to be quite small, i.e. 
~ 0.005 mag, consistent with our observational non-detection of 
any variability outside the eclipse. 

3.3.5. Spectral energy distribution 

We cross-identified SDSS1212-0123 with the database from 
the Galaxy Evoluiori Explorer (GALEX (iMartin et al.l 120051; 
iMorrissev et al.l l2005l) ). and found a detection in the far and 
near ultraviolet (FUV and NUV). The magnitudes are wipuv = 
16.79 + 0.03 mag and miMuv = 16.81 +0.02 mag, exposure times 
were 150 sec. FUV and NUV fluxes can provide an estimate of 
the effective temperature of the white dwarf for a certain log g, 
assuming that all the flux in the UV is emitted by the primary. 
White dwarf models for log g = 7.5 and log g = 8.0 and effective 
temperature in the range 6000 - 100000 K, were folded over the 
FUV and NUV filters. The calculated flux ratio FUV/NUV was 
compared with the observed for SDSS1212-0123. The GALEX 
flux ratio implies T^ff ~ 13000 K, significantly colder than what 
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Fig.9. Spectral energy disti-ibution of SDSS1212-0123. 
GALEX near and far ultraviolet and 2MASS infrared fluxes 
(black circles), optical SDSS spectrum (black line). A white 
dwarf model of T^ff = 17500 K and log^ - 7.5 (blue dashed 
line) and the spectum of LHS1504 with spectral type M5 from 
Legget's library (red dots) are shown for comparison . 



we obtain from the optical spectrum in Sec. 13.3.11 However, dis- 
crepant temperatures from GALEX UV and optical photometry 
were noticed earlie r from an analysis for a large number of white 
dwarfs (~ 250) bv lKawka & VennesI 120071) . This shows that one 
cannot expect the same UV and optical temperatures in a case- 
by-case basis, but at best on a statistical average. For the time 
being we accept the temperature from our fit to the SDSS spec- 
trum, which grossly reflects the UV to optical SED. 

The spectral energy distribution is shown in Fig.|9l including 
ultraviolet, optical and infrared fluxes from 2MA .SS. A model 
spect rum for a white dwarf of pure Hydrogen jKoester et alj 
|2005|) with effective temperature of 17500 K and logg = 7.5 
and a spectrum of the M5 star LHS 1504 from Legget's library^ 
are shown for comparison (iLeggett et al.ll2000h . 



Fig. 10. The ranges of masses of the white dwarf and the red 
dwarf coming from: the decomposition of the SDSS spectrum; 
the infrared brightness; the eclipse length for an inclination of 
90°, and, the detailed light curve modelling. Each of the areas 
is labeled correspondingly, the intersection of the four different 



methods occurs for M^f 



0.26-0.29 M„ andM, 



0.44-0.46 



http://ftp.jach.hawaii.edu/ukirt/skl/dM.spectra/ 



3.3.6. Binary parameters summary 

Figure [TO] shows the different ranges for the masses of the pri- 
mary and the secondary from the spectral decomposition fit 
(Sect. l3.3Tn i. the /T-band luminosity-mass relation (Sect. 13.3^ . 
the radial velocity amplitude and eclipse length (Sect. 13.33]) 
and the detailed light curve fitting (Sect. 13. 3^ . Of course, the 
different methods are not entirely independent, e.g. the con- 
straints from the eclipse length/radial velocities studies and the 
detail light curve fitting basically use the same information 
with the only difference being that we could derive a clear 
lower limit from the latter. The dark shaded region in Fig.fTOl 
represents the ranges of steUar masses in agreement with all 
the derived constraints i.e., M„d - 0.46 - 0.48 Mq, Msec = 
0.26 - 0.29 Mq, implying a radius of the secondary star in the 
range /^se c - 0.28 - 0.31 R© using the empirical M-R rela- 
tion from lBavless & Orosd 02006.) and R„i = 0.016 - 0.018 Rq 
(logg - 7.5 - 7.7). We adopt these values as the most probable 
ones and all finally accepted stellar and binary parameters based 
on Sloan-data, other catalogues and our own follow-up observa- 
tions are collected in Table [3] 



4. Evolutionary state 

The post CE evolution of compact binaries is driven by angu- 
lar momentum loss due to gravitational radiation and - perhaps 
much stronger - magnetic wind braking. Unfortunately, the lat- 
ter mechanism is currently far from being well constrained, and 
predicting and reconstructing the post CE evolution sensitively 
depends on the assumed prescription for magnetic braking. 

However, the di s rupted magnetic braking scenario proposed 
by iRappaport et al.l ( Il983h can still be considered the standard 
model for magnetic braking in close compact binaries. In this 
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Table 3. Stellar and binary parameters of SDSS1212-0123. 



Parameter 


Value 


Parameter 


Value 


R.A. (72000.0) 


12 12 58.25 


Teff (K) 


17700 ± 300 


Dec. (72000.0) 


-01 23 10.1 


log^ (dex) 


7.5 - 7.7 


u 


17.045 ±0.020 


Sp(2) 


M4± 1 


g 


16.769 ±0.013 


^wd (Ro) 


0.016-0.018 


r 


16.936 ±0.013 


^sec (Ro) 


0.28-0.31 


I 




M . CM "1 


4.6 - 4.S 


z 


16.136 ±0.018 


M,ec (Mo) 


0.26 - 0.29 


J 


14.83 ± 0.03 


Port (days) 


0.3358706(5) 


H 


14.35 ± 0.05 


K2 (kms-') 


181 ±3 




13.93 ±0.05 


72 (kms"') 


17 ±3 


mpuv 


16.79 ± 0.03 


a(Ro) 


1.8 ±0.1 


™NUV 


16.81 ± 0.02 


'mill 


82° 


d(pc) 


230 ± 20 







scenario it is assumed that magnetic braking ceases when the 
secondary star becomes fully convective at Msec ~ O.SM© 
(which corresponds to Porb ~ 3 hrs). Although observations 
of the spin down rates of single stars do drastically disagree 
with the predictions of disrupted magnetic braking dSills et all 
l200ft) . it remains the only consistent theory explaining the or- 
bital period gap i.e. the observed deficit of CVs in the range 
of Porb ~ 2 - 3hrs. Moreover, first results of our radial ve- 
locity survey of PC EBs seem to suppo rt the idea of disrupted 
magnetic braking ( ISchreiber et alj|2008l) . To predict and recon- 
struct the post CE evo l ution of SDSS 1212-0123 according to 
ISchreiber & Gansickd (l2003h . we therefore assume disrupted 

magnetic braking. 

First, we interpolate the cooling tracks of I Wood! (Il995h and 
estimate the cooUng age of SDSS1212-0123 to be 6.8 x 10^ yrs 
(see top panel of Fig.fTTTi. Second, according to the mass derived 
for the secondary star (Msec ~ 0.27Mo) we assume that, since 
SDSS 12 12-0123 left the CE phase, the only mechanism driving 
the evolution of SDSS1212-0123 towards shorter orbital peri- 
ods is (and has been) gravitational radiation. As shown in Fig.fTTI 
(bottom panel), SDSS1212-0123 left the CE phase with an or- 
bital period of Pce ~ 8.07 hrs, very similar to the present value. 
Significant changes in the orbital period are predicted to occur 
on timescales longer than the current cooling age of the white 
dwarf. In ~ 1.8 x 10'° years SDSS1212-0123 will eventually 
become a CV within the orbital period gap, however, giving that 
it's calculated PCEB lifetime exceeds the age of the Galaxy it is 
not representative of the progenitors of todays CV population. 

5. Conclusions 

From optical photometry we conclude SDSS1212-0123 is a 
eclipsing PCEB with an orbital period of 0.336 days and an 
eclipse length of 23 min. From spectroscopic follow-up obser- 
vations we have derived a systemic velocity of 17 + 3 km/s and 
a semi-amplitude of the radial velocity of 181 ± 3 km/s. From 
the SDSS spectrum we derived T^s = 17700 + 300 K, log^ = 
7.53 + 0.2 implying a mass in the range 0.33 - 0.48 and 
a secondary spectral type M4 + 1, and a distance to the system 
of 230 + 20 parsecs. From infrared photometry, using a mass- 
luminosity empirical relation we derived Msec - 0.26 + 0.03 M©. 
We have calculated the radius of the secondary star using an em- 
pirical mass-radius ralation. The mass function, combined with 
the eclipse length, points towards the high end of the allowed 
mass range of the primary, i.e. M„d ~ 0.46 - 0.48. We have 
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Fig. 11 . Top panel: Interpolating the cooling tracks from lWoodI 
(199$) and according to the current temperature of the white 
dwarf (Feff = 17700 K) we derive for SDSS1212-0123 a cool- 
ing age of ~ 7x10^ years. Bottom panel: Assuming gravitational 
radiation as the only angular momentum loss mechanism we re- 
constructed the post-CE evolution of SDSS1212-0123 and find 
that it left the CE phase with an orbital period of Pce ~ 8.07 hrs. 
Apparently, SDSS1212-0123 has passed only a small fraction 
of its PCEB lifetime and it will take ~ 1.8 x 10^'^ years until 
SDSS1212-0123 will become a CV. At that moment the white 
dwarf temperature will be Teff ~ 4000 K and the system will be 
inside the period gap (grey bar). 

modeled the / band light curve and find the inclination of the or- 
bit to be / > 82°, and the masses to be consistent with previously 
determined values. The different methods applied are all consis- 
tent with M„d = 0.46 - 0.48 Mq, implying /?wd = 0.016 - 0.018 
Rq (logg = 7.5 -7.7) for the primary and Msec = 0.26-0.29 Mq, 
/?sec = 0.28 - 0.31 R© for the secondary. We have reconstructed 
and predicted the post CE evolution of SDSS1212-0123, find- 
ing that SDSS1212-0123 at the end of the CE phase had a very 
similar orbital period. The only mechanism involved in shrink- 
ing the orbital period is and has been gravitational radiation. As 
the PCEB lifetime of SDSS1212-0123 exceeds the Hubble time 
we conclude that it is not representative of the progenitors of the 
current CV population. 
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